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> ■ ABSTRACT 

OO 



A theoretical formulation for the two-point correlation function on a light-cone is 



(N 

^^ ' developed in the redshift space. On the basis of the previous work by Yamamoto & 

Suto (1999), in which a theoretical formula for the two-point correlation function on 

, a light-cone has been developed in the real space, we extend it to the formula in the 

Q^' redshift space by taking the peculiar velocity of the sources into account. A simple 

Q . expression for the two-point correlation function is derived. We briefly discuss QSO 

^ , correlation functions on a light-cone adopting a simple model of the sources. 



Subject headings: cosmology: theory - dark matter - large-scale structure of universe 



1. Introduction 

The clustering of high-redshift objects is one of the current topics in the fields of observational 
cosmology and astrophysics. The high-redshift objects of z ^ 1 are becoming fairly common, and 
evidences of the clustering nature of such cosmic objects are reported in the various observational 
bands, e.g.. X-ray selected AGNs (Carrera et al. 1998), the FIRST survey (Cress et al. 1996; 
Magliocchetti et al. 1998), high-redshift galaxies (Steidel et al. 1998; Giavalisco et al. 1998), and 
QSO surveys (Groom &: Shanks 1996; Boyle et al. 1998). The statistics of these high-redshift 
objects are becoming higher, and we will be able to discuss the clustering at a quantitative level 
precisely in near future. From a theoretical point of view, the most important subject is to clarify 
the physical process of the formation history of these objects. The standard theoretical framework 
for the cosmic structure formation is based on the cold dark matter (GDM) model with gaussian 
initial density fluctuations. The clustering nature of the high-redshift objects provides us with 
many kinds of tests for the theoretical models (e.g.. Peacock 1998; Jing & Suto 1998). 

When analyzing the clustering nature of the high-redshift objects at a quantitative level, we 
must take the light-cone effect into account properly. Namely, such cosmological observations are 
feasible only on the light-cone hypersurface defined by the current observer. And the effect of the 
time-evolution of the sources, i.e., the luminosity function, the clustering amplitude, and the bias, 
contaminates an observational data. Thus this light-cone effect is especially important to discuss 
the three-dimensional two-point correlation function of the high-redshift objects. Some aspects 
of the light-cone effect has been discussed (Matarrese et al. 1997;Matsubara, Suto, & Szapudi 
1997; Nakamura, Matsubara, & Suto 1998;Laix & Starkman 1998). Recently one of the authors 
(K.Y.) & Suto developed a formulation for the two-point correlation function for the high-redshift 
objects defined on the light-cone hypersurface (Yamamoto & Suto 1999: hereafter Paper I). The 
expression for the two-point correlation function on the light-cone was derived in a rigorous 
manner starting from first principle corresponding to the conventional pair-count analysis. This 
investigation is very important because it gives a rigorous relation between an observational data 
processing and a theoretical prediction as to the two-point correlation function on a light-cone for 
the first time. However this investigation is restricted to the formula in the real space, though 
observational maps of the high-redshift objects are obtained in the redshift space. 

It is well known that the peculiar velocity of sources distorts their distribution in the redshift 
space (e.g., Davis & Peebles 1983; Kaiser 1987; Hamilton 1997). And this effect has been discussed 
as a probe of cosmological density parameters (e.g., Szalay, Matsubara, & Landy 1998; Nakamura, 
Matsubara, & Suto 1998; Matsubara & Suto 1996; Hamilton & Culhane 1996; Heavens & Taylor 
1995; Suto, et al. 1999). In the previous paper (Paper I) the effect of the redshift-space distortion 
due to the peculiar velocity of the sources is not taken into account because it was formulated 
in the real space. From a practical point of view, the formula in the redshift space must be 
developed. The purpose of the present paper is to develop such a theoretical formula for the 
two-point correlation function on a light-cone hypersurface by taking the redshift-space distortion 
due to the peculiar velocity into account. 



The paper is organized as follows: In §2 we develop a formulation for the two-point correlation 
function on the light-cone hypersurface in the redshift space in order to incorporate the linear 
redshift-space distortion. The expression for the two-point correlation function is presented in a 
rather simple form by using appropriate approximations. The main result is equation (32). As 



a demonstration of the usefulness of our formalism, we apply the formula to QSO correlation 
functions adopting a simple model of source distribution and cosmological models. The validity of 
the plane-parallel, or distant observer, approximation, for the correlation function of high-redshift 
objects is also discussed. §4 is devoted to discussion and conclusion. Throughout this paper we 
use the units in which the light velocity c is unity. 



2. Two-point correlation function in the redshift space 

In this section we develop a theoretical formulation for the two-point correlation function 
on a light-cone hypersurface in the redshift space by taking the peculiar motion of sources into 
account. In the present paper, we focus on the spatially-flat Friedmann-Lemaitre universe, whose 
line element is expressed in terms of the conformal time r] as 

ds^ = a\r,) [-dr,^ + dx^ + X^dflf^^] . (1) 

Here the scale factor is normalized to be unity at present, i.e., a(??o) = 1- The Friedmann equation 



''y = Hi(^ + a'nA, (2) 
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where Oa = 1 — Oq, the dot denotes ry differentiation, and Hq is the Hubble constant Hq = lOO/i 
km/s/Mpc. 

Since our fiducial observer is located at the origin of the coordinates (?? = %, X = 0); 
an object at x ^^^ V on the light-cone hypersurface of the observer satisfies a simple relation 
r] = r]Q — X- Then the (real-space) position of the source on the light-cone hypersurface is specified 
by (XjT)) where 7 is an unit directional vector. In order to avoid confusion, we introduce the 
radial coordinate r instead of X: ^^^ we denote the metric of the three-dimensional real space on 
which the observable sources are distributed, as follows, 

dslc = dr^ + r^dnf2)- (3) 

Denoting the comoving number density of observed objects at a conformal time r] and at a 
position (x,7) by n{r],x,7), then the corresponding number density projected onto the space (^) 
is obtained by 
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LC 



(r, 7) = n(?7, X, 7) 1,,^,,0-r, x~*r ■ (4) 

Introducing the mean observed (comoving) number density no(??) at time r] and the density 
fluctuation of luminous objects A(r/,x,7), we write 

n(,V,X,7) =no{r])[l + A{ri,x,7)], (5) 



then equation (^) is rewritten as 



nLC(r,7)=no(??)[l + A(7?,x,7)] \^^r,,-r, x^v (6) 

Note that the mean observed number density no(?7) is different from the mean number density of 
the objects rl{rj) at r; by a factor of the selection function (^(ry) which depends on the luminosity 
function of the objects and thus the magnitude-limit of the survey, for instance: no(??) = n{rj)4>{rj). 

In the similar way, if we know the peculiar velocity field, the corresponding quantity projected 
onto the space (y) is obtained by 

V^^{r, l) = Vciv, X, 7) \ri-.rio-r, x^r ' (7) 

where Vc{i],x,7) is the CDM velocity field. Here we assume that the peculiar velocity field of 
luminous objects agrees with the CDM velocity field. 

In Appendix A we summarized equations for the linear perturbation theory in the CDM 
dominated universe. Thus the linearized CDM density perturbation can be solved completely. 
However, the evolution of the source density fluctuations can not be solved completely since 
the bias mechanism is not well understood at present unfortunately. Then we must assume a 
model for the bias which connects the CDM density perturbations and the source number density 
fluctuations. In the present paper we assume the scale-dependent bias model: 



where Af^imirj) and 3^/^(1]) are the Fourier coefficients for the source number density fluctuation 
and the CDM density fluctuation, respectively (see also Appendix A). 

The next task is to describe the relation between the real space and the redshift space, since 
we consider the distribution of sources in the redshift space. First we consider how the peculiar 
velocity of a source distorts the estimation of the distance to the source. Let us assume that a 
source at redshift z (at a position (r, 7) in the real space) is moving with a peculiar velocity v. 
The observed photon frequency t'obs and the emitted photon frequency femit is related as 

z^obs = -^(l-7-^^)- (9) 

From this equation, we find the shift in the apparent redshift due to the peculiar velocity as 

6z = il + z){j-v) . (10) 

From the Friedmann equation (g), we have 

1 a^/^6z 



Combining equations ( |To| ) and (0), then we obtain the apparent shift in the comoving coordinate 
due to the pecuhar velocity as 

Sr = -5n = —^ ^ -v \ (12) 

where we defined 

^('^^ = /o 1 o i Y^ ■ ^^^) 

We introduce the variable s to denote the radial coordinate in the redshift space. Then a 
position in the redshift space is specified by (s,7), while the real space is done by (r, 7). The 
relation between the redshift position and the real position is 

s = r + 6r, (14) 

where 6r is specified by equation (p!2[). The conservation of the number of sources gives (Hamilton 
1997) 

ns(s,7)s dsdKl:^ = n {r,j)r drdO.^, (15) 

where ns{s,^) denotes the number density in the redshift space and n^^(r, 7) does the number 
density in the real space. These two equations (|l4[) and ([T5|) specify the relation between the 
redshift space and the real space. 

Now let us consider the two-point correlation function in the redshift space. We start from 
the following ensemble estimator for the two-point correlation function: 

Xs{R) = -jjLc / ~T^ / dsisldQ.;^-^ / ds2S2dn^^ns{si,^i)ns{s2,^2)S^^\si - S2 - R), (16) 

where si = (si,7i) and S2 = (82,72) and R = |R|, R = R/i?, and V^'^ is the comoving survey 
volume of the data catalogue: 

yLC ^ r-"- ^2^^ I ^^^ ^ ^{s^J - s^i„3), (17) 

with Smax = ^(zmax) and Smin = •s(zmm) being the boundaries of the survey volume. Equation 
dl^ ) is a natural extension of the ensemble estimator for the two-point correlation function in the 
redshift space (see also Paper I). 

By using equations (^) and ([I^), we rewrite equation {It) in terms of the variables in the 
real space: 

1 r dVi ~ r r 

'^s(^) = ^^ j -^ j dnrldVt;^^ j dr2rldVt;^,n^'^(ru^i)n^'^{r2,^2) 

x5(3)(xi+5xi-X2-5x2-R), (18) 



where xi + 5xi = (ri + 5ri,^i), X2 + 5x2 = (r2 + (Jr2, 72), and 6ri and 5r2 are given by (12). Then 
we approximate as 

d'-^^ (Xi + ,5X1 - X2 - 5X2 -R.)^(l + 5X1 • ^) (1 + 5X2 • T^) 5(3) (xi - X2 - R), (19) 
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where we can write 5x • d/dx = 6rd/dr since only the radial component of 5x has a non-zero 
value. By using this approximation and equation m) we derive the following equation from ([l8|): 
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yhC J ~^ 
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drirfdn^^ / dr2r^dn^^n^'^ {ri)n^'^ {r2) 



i=l 



d 

l + A(r„7i) + 5ri-- 
on 



5(3)(xi-X2-R), 



where we used the notations: 
n^^ir) = no(??) 
and Sri is understood as 

where i = 1, 2. 



V^Vo-r ' 



A(r,7) = A(r?,x,7) 1^- 



"70-'', X^'' ' 



ori = —5— 7 • ?; 



V^VO-rt, X^ri ' 



(20) 

(21) 
(22) 



Next we consider the ensemble average of the ensemble estimator Xs{R). Since 5ri and 5r2 
are the order of linear perturbation, then the ensemble average is written as 



{X,{R))=l({R)+WsiR), 



(23) 



where 



U{R) 



dn^ 
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^ dnrl /d% dr2rl dn^^n}^'^{n)n}^'^{r2)d^^\xi - X2 -K), 



(24) 



and 



m{R) 



dn. 



yLC J 47r 



^ / dnrf / d% / dr2rl / d^n^C^^^^^LC^^^) 



x( ( A(ri,7i) + Sn—^ (A{r2,j2) + Sr2-^ ) }S^^\xi - X2 - R) 



(25) 



In Appendix B we presented the explicit calculations for Ws{R). According to the result, 



equation (25) reduces to the following form within the linear theory of perturbation: 
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dridr2rir2nQ (ri)no (r2)£'i(r/o - ri)Di{rio - r2) 



V^^ttRj js 
2 
X / dkk^P(k) TT I A^i- -^ - -.^ - ^-2-0 I .•„ / 1.. /.2 ^ ^2 



^ 2 

J dkeP{k)Yl[bik; 7]o-r^)-k-^Vr] Jo[k^rf + r^2- 2nr2 cos 9) , (26) 



where 



^2 



d 



T^r = f{m - 0^ + fim ~ ^)l^ In ^r^no^(r) 1)1(7/0 - r)f{rjQ - r) 



d_ 
dr 



(27) 



and P{k) is the CDM power spectrum at present, Di{r]) is the hnear growth rate norniahzed to 
be unity at present, /(?]) is defined as /{rj) = dlnDi{rj)/dlna{i]), and cos9 should be replaced by 
cosO = (rf + rl — B?)/2rir2 after operating the differentiations with respect to ri and r2. 



Omitting the second term in the derivative (P7|), equation ([2q ) reduces to the simple form in 
the case R <C r^\a and R <^ r^nax (see Appendix B), 



Ws{R) 
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drr Uq (r) 



27r2 



1 + :j/3(r?o - r) + -P{m - rf 
6 5 



k^dkP{k) b{k-r]o-r)Di{rjo-r) 
MkR) , 



where (3{k;r]) is defined by 



/3(A:;77) 



fiv) 



1 dlnDi(?7) 



(28) 



(29) 



b{k;T]) b{k;7]) dlna{7]) 

and we assumed Smax = ''max (•Smin = ''minj- 

We can derive the following equation from a similar calculation in the above (see also Paper I): 

47r 



U{R) ~ 



r^drnoir). 



V^^ .n 

Following Paper I, we define the two-point correlation function on the light-cone hypersurface: 



(30) 



6^^(ii) 



{XsiR))-UiR) WsiR) 



U{R) 



U{R) 



(31) 



Substituting equations (Pq) and (30) into (|3l|), we have 



6^^(i?) 



27r2 



drr'^riQ^irf / drr^riQ^^r)'^ 
k'^dkP{k)b{k- r/0 - r)2L'i(r/o - r)^ 

jo(A;i?) . 



1 + o/3(^o - r) + -/3(??o - rf 
6 5 



(32) 



This is the final expression for the two-point correlation function on the light-cone hypersurface 
in which the linear redshift-space distortion is taken into account. Comparing this result with 
^^ (i?) in Paper I, the terms in proportion to I3{k; 7]) are the new terms which represent the effect 
of the linear redshift-space distortion. 



3. A Simple Demonstration 

In this section we apply the formula developed in the previous section to QSO two-point 
correlation function. Evidence for the spatial correlation in the QSO-distribution is reported 
(Croom &: Shanks 1996; Boyle et al. 1998), however, it seems difficult to draw definite cosmological 
conclusions from the comparison with the currently available data. Then we only demonstrate 
the usefulness of our formalism by calculating the QSO two-point correlation function based on a 
simplified model for the distribution and the bias model. As for the bias, we here consider the 
scale-independent bias model by Fry (1996): 

%) = l + 77^(feo-l), (33) 

where bo is a constant parameter. Note that the bias 6(77) at high-redshift becomes larger as 60 
becomes larger. Here we also assume that the sources are distributed in the range 0.3 < z < 3 
with a constant number density, i.e., uq = const. This model may be over-simplified, however, 
we have checked that the qualitative features have not been changed even when adopting more 
realistic models in Paper I. 

In figures 1 and 2 we show the two-point correlation function ^s (R) and other mass 
correlation functions for comparison. We show the case for the standard cold dark matter 
(SCDM) model in Fig. 1, in which we adopted Oq = 1, Oa = 0, /i = 0.5, and the CDM 
density power spectrum normalized as as = 0.56 (Kitayama & Suto 1997). The case for the 
cosmological model with a cosmological constant (ACDM model) is shown in Fig. 2, in which 
Oq = 0.3, Oa = 0.7, h = 0.7, o"8 = 1.0 are adopted. In each panel (a)-(c), the three lines show the 
correlation functions on a light-cone. The solid line shows ^s (R) of equation (p2|). The dashed 
line shows the case when neglecting the effect of the redshift-space distortion by setting /3 = in 
(I34 ). On the other hand, the panel (d) shows the linear and nonlinear mass two-point correlation 
functions defined on a constant time hypersurface z = 0, 1, and 2. 

From these figures it is apparent that the larger bias at the high-redshift derives the larger 
amplitude of the correlation function on a light-cone. Furthermore the effect of the redshift-space 
distortion always amplifies the correlation function from comparing the solid line and the dashed 
line in the panels (a)-(c), as expected. However the relative difference between the solid line and 
the dashed line becomes smaller as the bias becomes large and more effective. This is an expected 
feature because /?- factor becomes smaller as the bias becomes larger (see eg . ||29|| ) . 

We have also calculated the correlation function by adopting the exact expression (B20) 
instead of (28). The difference is less than 1 % for R ^ 100/i^^Mpc, and is negligible. Thus 
the formula (32) is a well approximated formula, and is an useful expession for the correlation 
function for high-redshift objects. The expression (32) is easily understood in an intuitive manner. 
Namely, the linear power spectrum in the redshift space is amplified by (1 + /3//k)^ over its 
unredshifted counterpart P{k) in the plane-parallel approximation, where /? is defined in the same 
way as (p9|) and /ik = 7 ■ k/Zc (e.g., Kaiser 1987; Hamilton 1997). This formula leads that the 



angle-averaged redshift power spectrum is amplified by the factor, (1 + 2/3/3 + /3^/5), over the 
unredshifted power spectrum. Thus (32) is the expected formula obtained by multiplying the 
factor, (1 + 2/3/3 + /?^/5), at each cosmological time over unredshifted counterpart P{k) in the 
correlation function in the real space. In this sense the formula (32) is based on the plane-parallel 
approximation.]^ And our investigation shows that the use of the plane-parallel approximation is 
valid for the (angle-averaged) correlation function of high-redshift objects on a light-cone. 

As we plotted absolute values of the two-point correlation functions in Figs. 1 and 2, then 
we can regard that the solid line and the dashed line show the anti-correlation at the large 
separation R ^ a few xlO/i^^Mpc in the SCDM model (see Fig. 1). Equation (32) implies that 
the zero-point of the correlation function is invariant even when the redshift space distortion is 
taken into account, as long as the bias does not depend on the scale k. The critical correlation 
length, where the correlation changes to the anti-correlation, is given by 

fi = iBM-'MP^ , ,34) 

r2o/iexp[— ilfe — V2hQo/^h] 

where we assumed the Harrison-Zeldovich initial density power spectrum and used the fitting 
formula for the transfer function (Bardeen et.al 1986; Sugiyama 1995). This critical correlation 
length may be observed in the upcoming 2dF and SDSS QSO surveys, and may be tested for the 
cosmological models and the theoretical models of bias. 



4. Summary and Discussion 

In this paper we have developed a theoretical formulation for the two-point correlation 
function for high-redshift objects on a light-cone in the redshift space. Our formula has been 
developed by extending the previous work (Paper I) to the formula in the redshift space. We 
have started our formulation from considering the ensemble estimator of the two-point correlation 
function in the redshift space, then we have calculated the ensemble average of the estimator. Thus 
our formula has been derived in a rigorous manner starting from first principle corresponding to 
the conventional pair-count analysis. The calculation was cumbersome, however, a rather simple 
expression (^21) has been derived. 



We have demonstrated the effect of the redshift-space distortion by showing the QSO 
two-point correlation function adopting a very simple model of the source distribution and the 
bias, though it seems premature to draw definite cosmological conclusions from comparison with 
currently available data. As discussed in the below, our model adopted in this paper may be 
oversimplified in order to compare with a real data sample. Nevertheless our investigation is 
instructive and we have shown how the redshift-space distortion affects the correlation function for 



^We thank T. Matsubara for his comment that our formula reproduces the formula inspired from such a 
consideration based on the plane-parallel approximation. 
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high-redshift objects on a light-cone (section 3). Our investigation shows that the redshift-space 
distortion becomes a smah effect for time-varying bias models which have large values at 
high-redshift. The validity of the plane-parallel approximation is also shown for the correlation 
function of high-redshift objects on a light-cone. 

There remain uncertainties to make precise theoretical predictions. First we did not attempt 
to examine possible bias models other than the model by Fry (1996). However, theoretical 
investigations for the time and scale dependent bias are just beginning (Fry 1996; Mo & White 
1996; Dekel &: Lahav 1998; Tegmark & Peebles 1998; Taruya, Koyama & Soda 1998). Conversely, 
the clustering of the high-redshift objects will be a good tool to test the bias models. Second 
we did not consider the realistic model for the time-evolution of number density in calculating 
QSO two-point correlation function. As for this point finite solution will be obtained in upcoming 
2dF and SDSS QSO surveys. Third we have only considered the linear theory of the density 
perturbations, and the nonlinear effect was not considered here. According to the previous work 
(Paper I), the non-linearity of the source density fluctuation becomes important only at small 
separation R ^ a, few /i^-'^Mpc (see also panel (d) in Figs. 1 and 2). And the effect of the 
non-linearity seems to be negligible at the large separation. However, the nonlinear effect in 
the redshift space has not been well understood especially for high-redshift objects, it must be 
investigated in future work. Probably numerical approaches would be needed for that purpose. 
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APPENDIX 



A. Review of the linear theory of the CDM density perturbations 

In this Appendix we summarize equations for the linear theory of the CDM density 
perturbations and explain the notations which are used in the present paper. The linearized CDM 
density perturbation in the CDM dominated universe obeys the following equations: 



6c + viii = 0, (Al) 

vi + -vi + ^1^ = 0, (A2) 



,-,« _L -v' -J- \T/N 

a 



v&i'|. = 47rGpa^5, = :^^<5c, (A3) 

where 6c is the CDM density contrast, vl is the CDM velocity field, ^ is the gravitational potential, 
which follows the gravitational poisson equation (^^), and \i denotes the covariant derivative on 
the three-dimensional space. 

As we are interested in the scalar perturbation, we expand the CDM density contrast 5c and 
the velocity field u* in terms of the scalar harmonics as follows (e.g., Kodama & Sasaki 1984): 

Sc{v,X,1) = r dkY,5tlMykUx.l), (A4) 

/■oo 

v"c{ri^X.l)= dkJ2vkimir])yii^{x,1), (A5) 

Jo , 

l,m 

where ykim is the normalized scalar harmonics: 

yMm{x,l) = Xi{x)YUn^), (A6) 

with 

Xiix) = \[^kji{kx), (A7) 

Yim{^^) and ji{x) are the spherical harmonics and the spherical Bessel function, respectively, k 
denotes the eigenvalue of the eigen-equation: ykim \i = —k'^ykim, and yii^ is defined as 

yiux,i) = -lykimix,i)^'- (A8) 



From the linearized perturbation equations (Al) — (A3), we have 



^kL + kvkim = 0, (A9) 

Vklm + -Vklm - k^klm = 0, (AlO) 

a 



S^qHq (c) 

2a '=''"' 



k'^kim = --^^^S',/^, (All) 
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where "^kim is the Fourier coefficient defined in the same way as ( [A 41) . Combining these equations, 
we have 

e + ^4t-|^C = 0. (A12) 

In the Friedmann-Lemaitre universe, the growing mode solution is well known: 

41(^)= 41(^0)^1 (a), (A13) 

with 

^■'°>=V!^+'-"°r''°'( ».+»-°(i-».) )'''- (A14) 

Here A is a constant to be determined so that Di is unity at present. From equations (|A9|) and 



(A5), we finally have 

vliri,X,l) = l^ dkJ2^-^^ykUx,l)^'. (A15) 



° Lm 



B. Calculation of Ws{R) 

In this Appendix we present an explicit calculation of Ws(i?): 

x^(^A(ri,7i) + 5ri^^(^A(r2,72)+<Jr2^)^5(3)(xi-X2-R). (Bl) 

Here A(r, 7) and (5r(r, 7) are explicitly written as 

A(r,7)=/ dkJ26'j^>^{rio)b{k;m-r)Di{rio-r)ykUrn), (B2) 

•^0 l,m 

, , ^. Z{r]o-r) 

Sr{r, 7) = v^ {r]o - r, r, 7) 

-no 

/•OO , , 

= / dfcE4/L(%)/(%-0^l(%-0^"'3^fc^n^(^,7)'^ (B3) 

•^0 l,m 

where we defined 

^ Zir^)dDr{v) 1 ^ ^ In 1^1(7?) 
•^^'^ Ho dri Di{r]) dlna{r]) ' ^ ' 

and we used equations @, (|T2|), (P), (|A13| ), and (|A15| ). Substituting equations ( p2|) and ( p3|) into 



equation (|B1|), we obtain 

1 /-df] 



^^^^^ = \kcJ^J ^^1^? / ^^71 / ^^2ri J dn^^ 
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72 V 



Zi ,m^ 



h,rn2 



d 



n b{h-7io-n)xl{n) + f{^o-n)k-^xl{np— ^(3) (xi - X2 - R 



i=l 



In addition, we use the relations 



5(3) (xi _ X2 - R) = — ^ /d^}^ g-*k.(xi-X2-R)^ 



-ikx 



(2vr) 

Z 771= — Z 



then equation (B5) becomes 

1 r dn 



Ws{R) 



R 



yLC J 47r 



driT]^ / dO^^ / dr2r2 / c^f^ 



''72 



xjdh E /^^2 E(4:u.(%)4:i:„..(%))^hm.(%)^:™.(^ 



72. 



h.mi 



l2,m2 



i=l 



.. d 



n Kh; m - ri)X'kM) + /(% - ri)ki'x'^^{ut' — 



dvi 



(2vr) 



^^Jd^in E (-i)^^iL,(A:ri)yL,Af,(0£)yz:,M,(^7i; 



"72. 



LiMi 

xAn E (^)''^JL2(fc^2)l2,Af,(^k)^i2Af2(^n 

1/2 Af 2 

x47r E ii)'''JLAkR)yL,Ahi^i.)YL,Msi^n) 
L3M3 



where A; = |k| and k = k/|k|. Integrating over i7^^, 0,^^ and il^ yields 

1 



>Vs(i?) 



yhC 



drir^ dr2r2nQ (ri)no (r2)£>i(??o - n)^i(??o - ^'2) 



Zi,mi Z2,m2 



i=l 



d 



n KA:.; % - n)x'^M) + /(770 - r.)A;-'xi;(r,)i'-» — 



du 



X l^js / rf'k(-z)'^-'^ J,, {kn)ji^ ikr2)JoikR)Yi^rn, i^i,)Yi;^, {n^: 



(B5) 

(B6) 
(B7) 



(B8) 



(B9) 
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And the further integration over Oj^ gives 



^«(^) = ^ J dnrf J dr2rln^'^{n)n^'^ir2)Di{7]o - n)D,{rjo 



r2 



X 



A;, 



d 



f dki /"dA:2^(2Z + l)P(A;i)J(A;i 
J J ^ 

X n Wk^; m - n)xl{n) + /(% - r,)kfxl{np-^ 

Xj^ I dkk'^ji{kri)ji{kr2)jo{kR) , 
where we used the relation of the gaussian random field in the linear theory: 

(4llimi(%)42*2m2(^0)) = P{kl)6{ki - k2)5l,l^6„^,m2, 

and 5i-^i^ and Smim2 ^^^ the Kronecker's delta. Integration by parts yields: 



(BIO) 



(Bll) 



X fdkiklP{ki)Y,{2l + l)^ f dkk^Ji{kn)ji{kr2)joikR). 

J ^ IT J 



i=l 



n ( Hh;vo - ri) - k^ ^Vr^ jjiikm) 



(B12) 



where we used (A7), and the operator T>r is defined by equation (|27|). Here the boundary terms 
are omitted since we can show that the boundary terms are the higher order terms of ©(iZ/rmax) 
or ©(iJ/rmin), which are negligible in the case R <C r^am and R ^ rmax- 

Noting the relation (e.g., Magnus et al. 1966): 



^ dkk^MkrMkr2)UkR) 



4f^^'(^%fe^) {\ri-r2\<R<n+r2), 
(i? < |ri -r2|, R> n +r2] 



(B13) 



we find 



Ws(i?) = :^7Lc^// dridr2rir2nl^{ri)nl^{r2)Di{7]o - ri)Di{7]o - r2] 



dk,klP{k,) ^(2/ + l)Pi ( "'+/' ^' 



n y>{ki]m-n) -k^'^VrAjiikiVi) 



(B14) 



where 5 denotes the region |ri — r2| < R <ri + r2. 
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By using the additional theorem for the spherical Bessel function: 

^(2/ + l)Pi{cose)ji{kri)ji{kr2) = JQ(k^ rl + rl- 2rir 2 cos d), 
I 

we can write 



(B15) 



m{R) 



1 1 

2 r 

dkk'^Pik) W 

i=l 



dridr2rir2nQ (ri)no (r2)-Di(??o - n)-Di(r?o - ^^2) 



h{k;rjQ - Ti) - k V,.^ 



Jo 



{k^rl + rl-2rir2Cose\ , (B16) 



where cosO is replaced by cos^ = (ri^ + r2^ — E?)/2rir2 after operating the differentiations by ri 
and r2. 



Introducing the notation z = \Jr\ + r^ — 2rir2 cos 0, we can show the formulas: 
k- 



drf 



Jo[kz) = r— (ri - r2 cos 9y - 



kz 



d'^ ■ n \ hikz) , „.2 Ji{kz 
-^jo{kz) = —^{r2-ricose) 



kz 



(B17) 
(B18) 



and 



92 g2 

k '^TT^TT^Joikz 
''2 



5^2 5^2Jo(A:zj- ^, 



(ri — r2COs6')^(r2 — ricos^)^ ^ ■; (^ri — r2COSl 



-4cos0(ri — r2 cos0)(r2 — ri cos 6) + (r2 — ri cos^)^ > 



.p(2™=..:,. 



(B19) 



Using these formulas and omitting the second term in the derivative, i.e., Vj- ~ f{7]o — r)d'^/dr'^, 
we have 



m{R) ^ 



1 1 

V^^J Js 

2 



dridr2rir2nQ (ri)no (r2) 



X / dkk^P{k) H [bik;rjo - rj)Di{r]o - Vj) 

X jo{kR) + P{k;r]Q - r2)IiR;ri,r2) + P{k;r]o - ri)I{R;r2,ri) 
+[3{k; 7/0 - ri)P{k; tjq - r2)J{R; n, r2) , 

where (3{k;r]) is defined by equation (p9|), and I{R;ri,r2) and J{R;ri,r2) are defined by 

/(i?;ri,r2) = ^-5 ^3— | } , (B21) 
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and 



J{R;ri,r2) 



h{kR) 



J4{kR) (R'^+rj-r^\'^( R'^ + r'^-rf^'^ 



{kRY 

k{k^ , , 
^ R^ \ 

h{kR) 



2<d±d^^\\, 



2rir2 



2ri 



2r2 



kR^ 



R'^ +r1- rl R^ + rl- rf rl + rl- R^^ 



n r2 2rir2 

respectively. 

Since we are generally interested in the case of i? ^ ?^max5 we can use the approximation: 



(B22) 



dridr2 ~ / dri I dx, 



(B23) 



where we introduced x = r2 — ri. By expanding I{R; ri, r2) and J{R; ri,r2) in terms of x, we have 

(B24) 
(B25) 



/(«;.,.,) = M^-A^.= 



kR 



i?2 



Integration by x leads to the final expression: 



,(ft,„,,.3MS_6M^,^ + «(^,. 
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Fig. 1. — Absolute values of the two-point correlation function for QSO on a light-cone and the 
linear and nonlinear mass two-point correlation functions in the standard CDM model, where 
we adopted Qq = 1, Q\ = 0, h = 0.5, and the CDM density power spectrum normalized 
as (Tg = 0.56 (Kitayama & Suto 1997). The parameter bo for the bias model are adopted as 
(a)6o = 1, (b)6o = 1-5, (c)6o = 2. Here we assumed that the sources are distributed in the range 
0.3 < z < 3 with a constant number density. In each panel (a)-(c) the solid line shows our ^s (R), 
the dashed line does the case when the redshift-space distortion is neglected by setting /3 = 0, (d) 
linear (lower curves) and nonlinear (upper curves) mass correlation functions by Peacock & Dodds 
(1996) defined on constant-time hypersurfaces z = 0, 1 and 2. 
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Fig. 2. — Same as Fig. 1 but for the ACDM model, in which we adopted 0,q = 0.3, Q\ = 0.7, h 
0.7, £78 = 1.0. 



